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PREFACE 
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Chief,  Structures  Laboratory  (SL),  and  James  T.  Ballard,  Assistant  Chief,  SL, 
and  Dr.  Jimmy  P.  Balsara,  Chief,  Structural  Mechanics  Division  (SMD),  SL. 

Dr.  Paul  F.  Mlakar,  formerly  of  SMD,  was  involved  in  the  planning  and  direct¬ 
ing  phases  of  the  work.  Instrumentation  support  was  provided  by  personnel  of 
the  Instrumentation  Services  Division,  WES.  Engineering  and  Construction 
Services  Division  personnel  fabricated  the  unusual  testing  system  used  to  load 
the  conduit  models.  Messrs.  Robert  E.  Walker,  C.  Dean  Norman,  and  R.  Stephen 
Wright,  SL,  provided  technical  guidance.  Mmes.  Linda  S.  Marble  and  Frances  M. 
Warren,  SL,  processed  the  test  data.  This  report  was  prepared  by  Dr.  Mlakar 
and  Mr.  Vincent  P.  Chiarito,  SMD.  The  report  was  edited  by  Ms.  Janean  C. 
Shirley,  Information  Products  Division,  Information  Technology  Laboratory, 

WES. 

The  previous  Director  of  WES  was  COL  Allen  F.  Grum,  USA;  the  present 
Commander  and  Director  is  COL  Dwayne  G.  Lee,  CE.  Dr.  Robert  W.  Whalin  is 
Technical  Director. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used 

(metric)  units  as  follows: 

Multiply 

in  this  report 

By 

can  be  converted  to  SI 

To  Obtain 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

feet 

0.3048 

metres 

inches 

2.54 

centimetres 

inch-pounds  (force) 

0.1129848 

metre-newtons 

kips  (force) 

4.448222 

kilonewtons 

kips  (force)  per  square  inch 

6.894757 

megapascals 

pounds  (force) 

4.448222 

newtons 

pounds  (force)  per  square  inch 

6.894757 

kilopascals 

square  feet 

0.09290304 

square  metres 

square  inches 

6.4516 

square  centimetres 

*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  the  following  formula:  C  =  (5/9)(F  -  32).  To  obtain  Kelvin  (K)  read¬ 
ings,  use  K  :  (5/9) (F  -  32)  +  273.15. 


STRENGTH  DESIGN  OF  REINFORCED  CONCRETE 
HYDRAULIC  STRUCTURES 

LOAD-MOMENT  CHARACTERISTICS  OF  REINFORCED 
CONCRETE  CIRCULAR  CONDUITS 

PART  I:  INTRODUCTION 

1.  The  investigation  summarized  herein  was  undertaken  to  determine  the 
effect  of  curvature  on  the  resistance  of  axially  flexurally  loaded  reinforced 
concrete  members  and  to  recommend  any  appropriate  modifications  to  normal 
strength  design  practice  for  circular  conduits  used  by  the  US  Army  Corps  of 
Engineers.  In  the  analysis  of  reinforced  concrete  (RC)  members  loaded  by 
combined  moment  and  thrust,  the  effects  of  initial  curvature  on  the  response 
of  RC  conduits  beyond  the  elastic  range  have  in  the  past  been  neglected.  The 
Corps  of  Engineers  constructs  many  curvilinear  conduits  with  curvatures 
considerably  greater  than  those  of  structural  members  encountered  in  usual 
practice. 

2.  A  survey  of  the  Corps'  existing  and  planned  conduits  was  conducted 

to  see  the  ranges  of  R/h  ,*  p  ,  p'  ,  d  ,  and  f  which  exist  for  these 

structures,  where:  R  =  radius  of  initial  curvature  to  middepth  of  section, 

h  =  overall  depth  of  section,  p  =  gross  tension  steel  ratio,  p1  =  gross 

8  8 

compression  steel  ratio,  d  =  effective  depth  of  tension  steel,  fy  =  yield 
strength  of  reinforcing  steel.  This  survey  is  shown  in  Table  1  and  sketches 
of  the  different  shapes  are  shown  in  Figure  1.  For  the  range  of  R/h  values 
listed,  Timoshenko's  curved  beam  theory  (Timoshenko  1941)  can  be  applied.  In 
Figure  2,  this  theory  indicates  that  maximum  stresses  computed  without  consid¬ 
ering  curvature  would  be  incorrect  (unconservative  by  at  least  10  percent)  if 
the  Corps  structures  behaved  elastically.  To  the  best  of  the  authors'  knowl¬ 
edge,  the  effect  of  curvature  on  an  inelastic  RC  beam  has  not  been  studied. 

3.  The  investigation  consisted  of  analytical  and  experimental  studies. 
The  analytical  study  involved  a  curved  beam  analysis  to  predict  capacities  of 
RC  sections  with  curvature  effects  as  well  as  an  application  of  a  nonlinear 
finite  element  (FE)  model  for  concrete.  The  experimental  study  involved 

*  For  convenience,  symbols  are  listed  and  defined  in  the  Notation 
(Appendix  E). 


For  Elastic  Materials  -  Error  Between  Linear 
and  Hyperbolic  Maximum  Stress  Value 
% 

35 

17 

10.9 

9.2 

3.2 


STRAIN  DISTRIBUTION 


Figure  2.  Comparison  of  stress  calculations, 
with  and  without  consideration  for  curvature 


one-dimensional  (1-D)  structural  and  two-dimensional  (2-D)  linear  elastic 
analyses  of  RC  conduit  models.  Three  RC  conduit  models  were  fabricated, 
instrumented,  and  tested.  The  results  of  the  experiments  are  plotted  as 
moment-thrust  values  and  are  compared  to  predicted  capacities  (interaction 
diagrams).  Results  are  discussed,  conclusions  are  drawn,  and  recommendations 
made . 

4.  This  investigation  provides  an  in-depth  look  at  the  flexural  re¬ 
sponse  of  circular  conduits,  but  applies  only  to  a  known  set  of  loading  condi¬ 
tions.  Although  the  response  of  conduits  is  a  soil-structure  interaction 
problem,  this  first  step  toward  its  solution  considers  only  the  design  load¬ 
ings  prescribed  in  EM  1110-2-2902  (Headquarters,  Department  of  the  Army  1969). 
Since  the  conduit  shapes  used  by  the  Corps  include  circular  and  noncircular 
geometries,  an  understanding  of  the  circular  structure's  behavior  is  a  logical 
prerequisite  before  analyzing  the  structural  response  of  more  complex  conduit 
shapes.  The  behavior  of  these  thick  elements  in  shear  must  be  considered  at 
some  future  time. 


p 

p 
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PART  II:  STRESS  AND  STRAIN  VARIATIONS  FOR  THE  STRENGTH 
DESIGN  OF  CURVED  CONCRETE  BEAMS 


5.  A  curved  beam  analysis  was  formulated  to  evaluate  the  response  of  RC 
conduits.  The  curved  beam  analysis  was  chosen  since  the  top  portion  of  each 
shape  for  a  conduit  represents  a  thickly  curved  beam.  There  is  very  good 
agreement  between  curved  beam  theory  and  the  elasticity  solution  (Airy  stress 
function)  (Timoshenko  1941)  for  stresses  in  a  thick  hollow  cylinder  under  a 
given  load  condition  as  shown  in  Figure  3-  The  Hognestad  constitutive  law  for 
concrete  is  used  to  investigate  the  inelastic  response  of  thick,  RC  curved 
beams.  The  curved  beam  analysis  shows  that  the  present  American  Concrete  In¬ 
stitute  (ACI)  practice  for  determining  flexural  resistance  of  RC  straight  beam 
members  (ACI  1983)  is  appropriate  for  the  Corps'  RC  curved  members. 

Curved  Beam  Analysis 


6.  Consider  the  curved  beam  subjected  to  a  positive  bending  moment  and 
axial  thrust  in  Figure  4.  Notation  for  Figure  4  is  as  follows: 
b  =  width  of  section 
c  =  depth  of  section  in  compression 

r  =  radius  of  initial  curvature  to  neutral  axis  of  section 
d<t>  =  small  angle  between  two  neighboring  cross  sections 
Ad4>  =  small  angle  of  rotation  due  to  moment  and  thrust 
The  distance  from  the  neutral  axis,  y  ,  is  positive  in  the  compressive  region 
and  negative  in  the  tensile  region.  If  it  is  assumed  that  plane  sections  of 
the  beam  initially  remain  plane  under  the  action  of  combined  moment  and 
thrust,  it  then  follows  that  the  strain  distribution  when  the  inner  surface  is 
in  tension  i3  given  by  (Timoshenko  1941): 

e(y)  .  y  .  a**  (1) 

r  +  y  d$  '  u 

where  e  =  strain.  In  this  hyperbolic  strain  distribution,  compressive 
strains  are  positive  values.  Similarly,  when  the  outer  surface  is  in  tension 
(subjected  to  a  negative  moment)  the  strain  distribution  is: 


e(y) 


y  Ad4> 
r  -  y  d«t> 


(2) 


where  positive  y  is  now  measured  from  the  neutral  axis  (N.A.)  toward  the 
center  of  curvature.  The  term  Ad<t>/d$  in  Equations  1  and  2  is  a  constant 
that  is  evaluated  by  setting  e(y)  equal  to  eu  at  y  equal  to  c  . 


Therefore,  Ad<t>/d4>  equals  e 


where  the  plus  sign  is  used  for  a 


(r  t  c) 
u  c  ’ 

positive  moment  section  and  the  minus  sign  is  used  for  a  negative  moment 
section.  It  is  further  supposed  that  the  stress-strain  relationships  for 
concrete  in  compression  (Hognestad,  Hanson,  and  McHenry  1955)  and  for  steel 
reinforcement  in  tension  and  compression  are  as  shown  in  Figure  5,  where: 

Eg  =  modulus  of  elasticity  of  steel  (29,000,000  psi*) 
fQ  =  concrete  stress 

f^  =  28-day  compressive  strength  of  standard  6-  by  12-in.  concrete 
cylinder 

f”  =  compressive  strength  of  concrete  in  reinforced  concrete  members 
(0.85  f'c) 

f  =  steel  stress;  for  this  study  GR60  (grade  60  steel)  is  used 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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e  =  strain  at  which  the  maximum  concrete  stress  f"  occurs  (0.002) 
o  c 

eu  =  limiting  useful  concrete  strain  (0.003) 

Finally,  in  accordance  with  accepted  procedure,  assume  the  tensile  strength  of 
concrete  to  be  negligible. 


b.  Steel  in  tension  or  compression 

Figure  5.  Stress-strain  relationships 

7.  Under  these  assumptions  the  load-moment  characteristics  of  a  given 
initially  curved  reinforced  concrete  cross  section  can  be  obtained  as  shown  in 
Figures  6  and  7  through  the  following  steps  adopted  from  an  analysis  for  an 
initially  straight  member  (Pfrang,  Siess,  and  Sozen  1964).  For  Figures  5  and 
7  Cs  ,  Cc  ,  and  Ts  are  defined  below: 

Cg  =  Resultant  force  in  the  compressive  steel 
Cc  =  Resultant  compressive  concrete  force 
Tg  =  Resultant  force  in  the  tensile  steel 

M  =  Moment  applied  to  section 

( "+"  indicates  positive  moment,  negative  moment 


N.A.  =  Neutral  Axis 
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a.  Define  the  shape  of  a  hyperbolic  strain  distribution  over  the 
depths  of  the  section  by  assuming  a  neutral  axis  location  r  in 
Equation  1  for  positive  moment  or  Equation  2  for  negative 
moment. 

b.  Determine  the  concrete  stress  distribution  consistent  with  the 
assumed  strain  distribution  through  the  constitutive  relation  of 
Figure  2a  and  calculate  its  resultant  axial  thrust. 


b  /  fc(E,y)  dy  = 


and  moment  about  the  neutral  axis 


b  f  yf„  (e,y)  dy  =  M 


where  fc(e,y)  is  given  by  the  relationship  illustrated  in 
Figure  5a.  The  governing  equations  for  Figure  5a  are 


fo(*«y)  =  fc 


for  o  <  e  <  e„ 


2e(V)  _  |eMJ2 


fc(«*y)  ■  fc 


for  eQ  <  e  <  eu  and  e(y)  defined  by  Equations  1  and  2. 


c.  Determine  the  reinforcement  forces,  Cs  and  Ts  ,  correspond¬ 
ing  to  the  assumed  strain  distribution  by  the  constitutive  rela¬ 
tion  of  Figure  5. 

d.  Calculate  the  centric  thrust,  P  ,  and  the  bending  moment,  M  , 
acting  on  the  cross  section  which  are  statically  equivalent  to 
this  distribution. 


8.  If  a  series  of  such  strain  distributions  is  established  whose  magni¬ 
tude  corresponds  to  a  failure  of  the  cross  section,  then  the  resulting  (M,  P) 
points  are  the  ultimate  interaction  diagram.  For  purposes  of  practical 
design,  a  failure  occurs  in  tension  if  the  tensile  reinforcement  strain  es 
equals  its  tensile  yield  strain  -fy/eg  ;  or  in  compression,  if  the  maximum 
concrete  strain  ec  equals  the  ultimate  concrete  strain  eu  and  the  middepth 


concrete  strain  is  less  than  eq  ;  or  in  compression,  if  the  middepth  concrete 

strain  equals  e  and  is  less  than  e  . 

^  o  c  u 
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Results 


9.  The  resultants  Cc  and  Mc  of  the  concrete  stress  distribution 


in  an  initially  curved  cross  section  can  be  determined  by  closed-form  integra¬ 
tions  under  the  postulates  and  procedures  of  the  preceding  section. 

10.  Consider  first  a  member  loaded  by  a  positive  moment  which  causes 
compression  on  the  cross  section's  outer  surface  as  shown  in  Figure  6.  For  a 
failure  governed  by  the  yield  of  the  tensile  steel 
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where  x'  =  c  -  h  .  For  a  failure  governed  by  the  middepth  strain  equal  to 
eQ  the  resultants  are  given  by  Equations  7  and  8  with: 
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and  substituting  k'  for  k  . 

11.  A  member  subjected  to  a  negative  moment  behaves  as  in  Figure  7. 
Cc  and  Mc  are  given  by 
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For  a  failure  in  which  the  middepth  strain  equals  eq  ,  Cq  and  Mc  are 
found  by  substituting  k'  from  Equation  9  for  k^  in  Equations  12  and  13. 


Discussion  of  Analysis 


12.  To  assess  the  implication  of  the  above  results,  one  must  know  the 
range  of  curvatures  and  other  variables  encountered  in  Corps  conduit  construc¬ 
tion.  To  this  end,  recent  Corps  publications  revealing  variable  values  for 
conduit  structures  were  quickly  surveyed  as  summarized  in  Table  1.  If  this 
survey  is  taken  to  be  representative  of  Corps  practices,  it  would  seem  that 
the  values  of  dimensionless  curvatures  R/h  range  from  2.5  to  3.8. 

13.  The  effect  of  the  most  severe  curvature  in  this  domain  is  super¬ 
imposed  on  a  previously  published  (Mattock,  Kriz,  and  Hognestad  1961)  analysis 
of  concrete  stress  distribution  in  Figure  8.  In  this  figure,  a  dimensionless 
measure  of  the  concrete  resultant  thrust  is  given  by: 
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Figure  8.  Properties  of  concrete  stress  distribution  at 
ultimate  strength  determined  from  tests  on  plain 

concrete  (Mattock,  Kriz,  and  Hognestad  1961) 


The  equations  describing  the  solid  curves  in  Figure  8  are  presented  in  Appen¬ 
dix  D.  The  dashed  curves  using  Equations  5  through  13  are  not  a  function  of 
the  cylinder  strength.  There  is  good  agreement  between  the  analytical  results 
of  this  study  for  R/H  ♦  =  and  independently  acquired  experimental  results 
for  members  with  no  initial  curvature.  It  is  also  noted  that  the  effect  on 
k  due  to  the  most  severe  initial  curvature  in  Corps  practice  is  within  the 
scatter  of  carefully  conducted  laboratory  experiments  on  uncurved  members. 

The  solid  lines  suggest  that  the  present  ACI  S1  factor,  which  was  shown  as 
the  "Proposed  Value"  in  Figure  8,  is  about  as  appropriate  for  Corps  practice 
as  for  normal  ACI  practice.  The  dashed  curves  for  the  most  extreme  case  of 
curvature  at  the  springing  line  section  do  lie  slightly  below  the  scatter 
between  3,000  and  5,000  psi  cylinder  strengths.  Results  from  the  tests  on 


the  model  conduits  will  be  used  to  evaluate  further  using  the  ACI  factor 

for  Corps  practice  for  designing  RC  circular  conduits. 

14.  Determining  cross-section  resistance  of  straight  members  under  com¬ 
bined  moment  and  thrust  is  compared  to  the  capacity  of  Corps  curved  members. 
This  can  be  seen  from  the  interaction  diagrams  of  Figures  9  and  10  for  cross 
sections  representative  of  Table  1.  In  these  figures,  the  resistance  of  a 
member  initially  curved  at  R/h  =  2.5  is  shown  by  the  symbols.  It  is  seen 
that  the  hyperbolic  strain  distribution  induced  by  the  initial  curvature 
causes  a  smaller  resistance  in  negative  moment  than  in  positive  moment.  How¬ 
ever,  when  these  resistances  are  compared  to  the  resistance  of  a  straight  beam 
(R/H  -*  °°  corresponds  to  a  straight  beam),  a  decrease  is  only  seen  near  the 
balance  point,  and  this  difference  is  within  10  percent. 


Figure  9.  Interaction  diagram,  symmetrically  reinforced 
cross  section  representative  of  Table  1 


15.  The  nonlinear  FE  program  NONSAP  was  used  to  study  the  behavior  of 
a  plain  concrete  curved  beam.  An  elastic-plastic  hardening  model  which 
accounts  for  fracturing  in  the  nonlinear  material  it  represents  and  assumes 
small  strains  in  rotation  was  used  to  represent  the  plain  concrete  (Chen  and 
Chen  1975).  It  should  be  noted  that  this  model  only  represents  the  portion 
of  the  concrete  stress-strain  curve  up  the  maximum  stress,  f"  .  In  Fig¬ 
ure  11,  the  Chen  and  Chen  model  and  the  Hognestad  stress-strain  curve  are 
compared.  The  symbols  depict  results  computed  from  the  application  of  the 
Chen  and  Chen  model  to  a  simply  supported  beam  under  combined  moment  and 
thrust  loading.  The  stress-strain  curves  in  Figure  11  are  both  for  uniaxial 
stress  conditions. 

16.  The  differences  noted  between  the  two  curves  in  Figure  11  are  the 
maximum  concrete  stresses  allowed  and  the  strain  value  at  which  the  maximum 
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Figure  11.  Comparison  of  the  Hognestad  stress-strain  curve  (solid) 
with  uniaxial  nonlinear  FE  data  (Chen  and  Chen  model,  depicted 

by  symbols) 

stress  occurs.  In  the  Hognestad  model,  there  is  a  factor  that  controls  the 
value  of  the  maximum  stress. 

17.  A  nonlinear  FE  analysis  was  performed  on  a  plain  concrete  curved 
beam  to  verify  the  validity  of  a  hyperbolic  strain  distribution  for  loading  in 
the  inelastic  range.  In  Figure  12,  the  grid  is  shown  for  this  analysis.  The 
grid  consisted  of  192  total  elements,  of  which  112  were  2-D  linear  elements  in 
the  outer  thirds  and  80  were  2-D  nonlinear  elements  located  in  the  middle 
third  of  the  plain  concrete  member.  The  curved  beam  was  simply  supported  at 
the  geometric  centroid  of  the  section.  Horizontal  load  couples  placed  at  the 
ends  of  the  beam  applied  a  constant  moment,  and  an  external  pressure  applied 
on  the  extrados  of  the  curved  beam  produced  a  constant  thrust  throughout  the 
beam.  The  two  load  curves  (force  couple  and  pressure)  were  defined  so  that 
the  neutral  axis  remained  close  to  the  edge  such  that  the  entire  cross  section 
was  in  compression.  This  simulated  the  procedure  used  by  the  Portland  Cement 
Association  (PCA)  to  conduct  eccentric  column  tests  to  determine  stress-strain 
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Figure  12.  Grid  used  for  the  NONSAP  analysis  of  RC  curved  beam 

(undeformed  grid) 


characteristics  of  reinforced  and  plain  concrete  rectangular  sections 
(Kaar,  Hanson,  and  Capell  1977;  Kaar  et  al.  1978).  In  Figure  12  the  loading 
and  boundary  conditions  used  in  this  analysis  are  illustrated. 

18.  Figure  13  shows  the  verification  of  the  hyperbolic  strain  dis¬ 
tribution  by  the  nonlinear  FE  model,  up  to  a  strain  of  0.0015  in. /in.  com¬ 
pression.  The  strain  profiles  are  from  the  middle  vertical  cut  of  the  beam. 
The  curve  marked  by  t  s  20  (load  step  20)  depicts  elastic  behavior  of  the 
section  whereas  the  curve  marked  by  t  =  50  depicts  a  portion  of  the  section 
behaving  plastically  near  the  edge  of  the  section.  At  t  =  20  ,  the  pressure 
is  495  psi  and  the  horizontal  load  couple  is  (4.424  in.)  (0.85  *  10"*  lb) 

=  376,890  in. -lb.  At  t  =  50  ,  the  pressure  is  750  psi  and  the  horizontal 
load  couple  is  (4.434  in.)  (0.114  x  10^  lb)  =  505,476  in. -lb. 

19.  Also  listed  in  Figure  13  are  the  internal  forces  denoted  with 
subscript  I  and  the  applied  forces  denoted  with  subscript  A.  The  internal 
forces  were  calculated  from  the  stress  output  and  the  applied  forces  were  the 
horizontal  load  couple  and  the  thrust  due  to  uniform  pressure.  The  comparison 
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between  the  internal  and  applied  forces  shows  good  agreement.  Also  shown  is  a 
comparison  between  the  hyperbolic  and  linear  fits  to  the  numerical  results. 

p 

The  correlation  coefficient  squared,  rc  ,  is  listed  for  each  load  level. 
(Refer  to  Appendix  C  for  the  statistical  analysis.) 
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PART  III:  STRUCTURAL  ANALYSIS  OF  MODEL  RC 
CIRCULAR  CONDUITS 


20.  A  1-D  structural  analysis  was  performed  to  compute  the  internal 
forces  of  a  one-quarter  symmetrically  loaded  circular  conduit.  Three  loading 
cases  were  analyzed:  continuous  normal  pressure  (Anderson,  Haelsig,  and 
Reifel  1966),  the  3:1  Engineer  Manual  (EM)  design  load  condition  (Head¬ 
quarters,  Department  of  the  Army  1969),  and  a  five-point  load  simulating  the 
normal  load  component  of  the  3:1  EM  design  load  condition.  The  3:1  EM  loading 
is  shown  in  Figure  14.  The  equations  for  the  internal  moments,  thrusts,  and 
shears  are  shown  in  Appendix  A. 

21.  A  closed-form  solution, 

Airy  stress,  of  the  stresses  with  a 
given  wall  thickness  was  published 
in  an  Air  Force  Weapons  Laboratory 
(AFWL)  report  (Anderson,  Haelsig, 
and  Reifel  1966).  The  equations  for 
the  resultant  moments,  thrusts,  and 
shears  from  the  stress  distributions 
due  to  the  continuous  normal  pres¬ 
sure  distribution  are  also  presented 
in  Appendix  A. 

22.  In  Figure  15,  the  2-D 
conduit  is  idealized  by  the  1-D 
linear  elastic  line  structure.  The 
one-quarter  symmetric  structure  is 
1  degree  indeterminate;  therefore, 
both  the  equilibrium  equations  and  a  geometric  compatibility  condition  were 
required  to  solve  for  the  internal  forces.  All  of  the  material  properties 
were  assumed  to  be  represented  by  the  line  structure  at  the  geometric  center 
of  the  section.  (Actually,  it  is  more  accurate  to  assume  that  these  prop¬ 
erties  are  concentrated  at  the  neutral  axis,  but  the  neutral  axis  is  not  known 
a  priori.)  Only  the  normal  loading  can  be  modelled  experimentally  in  this 
study.  If  an  actual  soil  is  used,  the  3:1  EM  loading  can  be  more  accurately 
modelled  such  that  shear  transfer  exists  at  the  conduit-soil  interface. 


Figure  14.  EM  3:1  loading 
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1-D  idealization  of  RC  conduit 


23.  Comparisons  between  the  EM  and  the  two  AFWL  loadings  (1-D  and 
Airy)  are  presented  in  Figure  16.  The  analyses  were  normalized  with  respect 
to  the  applied  maximum  pressure  and  the  radius  of  curvature  to  middepth.  The 
2-D  results  were  obtained  from  a  problem  with  R  =  15  in.  and  h  =  6  in. 
Comparing  the  EM  loading  with  the  AFWL  1-D  loading,  the  EM  loading  gives 
50  percent  greater  moment  and  shear,  29  percent  more  thrust  at  the  springing 
line,  and  40  percent  less  thrust  at  the  crown.  Comparing  the  Airy  stress 
solution  with  the  AFWL  1-D  solution,  the  largest  differences  between  the 
thrusts,  shears,  and  moments  are,  respectively,  20,  20,  and  100  percent 
(at  0  =  45  deg,  where  the  1-D  analysis  gives  zero  moment,  and  the  Airy  stress 
solution  gives  -0.009  wR^  for  moment). 


PART  IV:  EXPERIMENTAL  STUDY  OF  THICK-WALLED 
CIRCULAR  CONDUITS 


24.  An  experimental  study  was  performed  to  verify  a  conventional  1-D 
structural  analysis  of  thick-walled  RC  circular  conduits.  The  curved  beam 
analysis  described  in  Part  II  was  used  to  predict  the  ultimate  capacity 
curves  of  the  curved  beam  RC  sections.  In  this  part  of  the  report,  only  the 
experimental  results  will  be  presented;  the  test  results  are  discussed  in 
Part  V.  Failure  loads  were  predicted  for  each  of  the  experiments. 

25.  In  addition  to  the  experimental  study  a  structural  analysis  of  a 
circular  conduit  subjected  to  the  EM  design  loading  (Headquarters,  Department 
of  the  Army  1969),  a  continuous  pressure  loading,  and  a  multipoint  loading  was 
performed.  The  pressure  and  the  test  loading  simulate  normal  components  of 
the  3:1  EM  design  loading  condition.  The  multipoint  loading  represents  loads 
by  hydraulic  rams  positioned  radially  about  a  circular  conduit  test  specimen. 


Test  Descriptions 


26.  The  purpose  of  the  tests  in  the  experimental  study  was  to  evaluate 
the  behavior  and  failure  modes  of  RC  circular  conduit  with  different  initial 
curvatures.  Three  RC  conduit  models  were  fabricated.  Two  were  identical 
(with  the  same  section  design  and  curvature  (R/h  =  2.5))  and  the  third  dif¬ 
fered  only  in  curvature  (R/h  =  4.0).  The  test  specimens  represent  a  scaled 
slice  of  a  circular  RC  conduit.  Cross  sections  of  Cl,  C2,  and  C3  are  shown  in 
Figures  17  and  18.  The  circular  shape  was  chosen  since  it  is  the  least  com¬ 
plex  and  the  tests  will  not  be  affected  by  any  abrupt  change  in  geometry  that 
might  complicate  interpreting  the  observed  behavior.  Also  the  first  series  of 
tests  establishes  a  standard  testing  procedure  for  future  tests  of  other 
conduit  shapes. 

27.  The  section  designs  of  the  models  are  representative  of  the  initial 
curvatures,  steel  percentages,  and  materials  of  the  range  of  prototype  struc¬ 
tures  surveyed,  as  listed  in  Table  2. 

28.  Normal-weight  concrete  was  used  with  a  maximum  size  aggregate  of 
1/2-in.  sand  and  normal  type  cement.  Steel  reinforcement  was  modelled  using 
No.  2  and  No.  3  Grade-60  rebars.  A  symmetrically  reinforced  section  was 
designed  for  the  tests. 
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load  channels  were  connected  to  50,000-lb  Precision  Universal  Interface 
ad  cells,  one  mounted  on  each  rod  end  of  the  hydraulic  jacks.  The  load 
11s  measured  the  actual  load  distribution  applied.  All  data  were  recorded 
two  32-track  magnetic  tapes.  Gage  locations  are  shown  in  Figure  19. 
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Figure  21.  Plan  view  of  testing  apparatus 


PART  V:  TEST  RESULTS 


Material  Properties 

32.  For  analyzing  the  test  results,  the  following  material  properties 
were  required:  28-day  uniaxial  compressive  concrete  strength;  uniaxial  com¬ 
pressive  concrete  strength  at  day  of  test;  and  yield,  ultimate,  and  rupture 
value  of  the  steel  reinforcement.  In  Table  4  the  results  from  the  uniaxial 
compressive  tests  and  the  tensile  tests  of  the  steel  reinforcement  are  listed. 
The  design  material  properties  representative  of  values  from  Table  1  for  the 
concrete  28-day  strength,  f^  ,  and  the  steel  reinforcement  yield  strength  are 
4  and  60  ksi,  respectively  (remembering  that  new  common  practice  is  to  use 
Grade-60  steel  instead  of  Grade  40).  Three  standard  6-  by  12-in.  concrete 
cylinders  were  tested  at  28  days.  Very  good  control  was  achieved  for  a  three- 
cylinder  average  28-day  strength  of  3,980  psi.  The  other  three  dates  in 
Table  4  were  the  days  that  models  Cl,  C2,  and  C3  were  each  first  tested. 

33.  The  values  of  f^  on  the  test  day  are  used  as  parameters  for 
the  concrete  material  property  according  to  Hognestad's  constitutive  law 
(Hognestad,  Hanson,  and  McHenry  1955).  In  models  Cl  and  C2,  a  weighted  aver¬ 
age  was  performed  to  compute  an  appropriate  value  for  the  yield  strength, 

fy  ,  of  the  steel  reinforcement. 

34.  Test  cylinders  and  models  were  cured  in  separate  buildings  because 
of  the  convenience  of  the  location  of  appropriate  test  equipment.  However,  it 

O 

is  noted  that  the  average  temperature  difference  probably  ranged  from  5  to 
10°  F  (i.e.,  the  temperatures  would  be  80°  to  85°  F  compared  with  85°  to 

O 

95  F).  It  is  not  known  exactly  how  much  of  a  difference  this  makes  in  the 

actual  f^  for  the  model.  However,  there  are  effects  other  than  temperature 
that  account  for  actual  differences  between  f^  for  test  cylinders  and  the 
models  (ACI  1970).  It  is  assumed  in  tnis  study  that  this  difference  is  neg¬ 
ligible  compared  to  other  experimental  errors  encountered. 

Test  Procedure 

35.  Before  the  first  specimen  was  tested,  a  load  test  was  conducted  to 
check  the  prescribed  load  distribution.  A  very  stiff  circular  specimen  with 
the  same  outside  dimensions  as  the  model  conduits  was  tested  with  the  same 


setup  proposed  for  the  conduit  tests.  The  results  from  the  load  verification 
are  presented  in  Table  5  and  agree  closely  with  the  predicted  loads.  The  mea¬ 
sured  loads  provided  a  check  on  the  load  distribution. 

36.  After  the  load  distribution  was  set,  a  model  conduit  was  placed  in 
the  reaction  frame  such  that  the  center  line  of  the  ram  aligned  with  the 
center  line  of  the  model.  The  model  was  supported  at  four  points  by  steel 
plate  shims  and  two  4  by  4  by  1/2-in.  pieces  of  Teflon.  The  Teflon  pads  were 
cleaned  with  a  solvent  beforehand  and  provided  a  very  good  low-friction  sur¬ 
face  to  allow  free  translation  of  the  model  during  loading.  With  the  model  in 
place,  the  model  instrumentation  channels  were  connected. 

37.  Before  loading,  each  rod  was  positioned  to  be  just  touching  load- 
bearing  plates  positioned  against  the  model.  The  pressures  in  the  jacks  were 
simultaneously  increased  and  all  the  data  channels  were  recorded  on  two 
32-track  magnetic  tape  machines.  The  loading  rate  for  each  of  the  tests  was 
approximately  100  psi/min  (2,000  lb/min  in  the  jack  with  the  largest  load). 
Cracking  was  detected  by  ear,  but  was  not  easily  seen.  Immediate  output  was 
in  the  form  of  an  x-y  plot  (one  load  channel  used  as  a  reference  versus  any 
two-strain  channels)  and  strip  charts  for  all  channels.  The  strip  charts  gave 
an  approximate  look  at  numerical  results  and  indicated  which  gage  channels 
were  ineffective  (either  because  of  noise  or  a  bad  gage).  The  data  stored  on 
the  32-track  tapes  were  then  digitized  and  placed  in  computer  files  accessible 
through  time-sharing  terminals. 

Observations 

38.  The  first  two  specimens,  Cl  and  C2,  remained  virtually  intact  when 
tested  with  the  full  capacity  of  the  loading  system.  Cracking  was  detected  by 
ear  during  the  tests  but  no  cracks  or  distress  in  Cl  or  C2  could  be  seen  with¬ 
out  closely  inspecting  each  after  the  tests  were  complete.  Retests  on  both 
specimens  showed  similar  results;  however,  more  tensile  cracks  were  visible. 
Hence,  Cl  and  C2  were  not  tested  to  failure  under  the  prescribed  load  distri¬ 
bution.  The  first  yield  was  expected  to  occur  at  the  crown  and  the  second  at 
the  springing  line  for  a  dominant  flexural  failure  mode.  The  crack  patterns 
observed  are  shown  in  Figure  22. 

39.  Model  C3  was  tested  to  failure  and  the  progression  of  the  collapse 
mechanism  being  formed  was  observed.  First  yield  and  plastic  hinge  were 


\  --  v  ■ 


h  =6" 
R/h  =  2.5 


.  SOME  SIGNS  OF  COMPRESSION 
CONCRETE  BEGINNING  TO  SPALL 


FROM  2  PT  LOADING  OF  SPECIMEN 
Cl  FOR  CALIBRATION  CHECK  AFTER 
1ST  TEST  ON  SPECIMEN  Cl 


- MINOR 
CRACKS 
ON  SURFACE 


■  INCIPIENT  CRUSHING 
OF  CONCRETE  VISIBLE 


-NO  VISIBLE 
CRACKS 


•  4 

■-  •*.  •*.  ■*.  •% 
f  v  *.«tv  V, 

V. 

AV.V.V.V. 


''-.•  V  W' 

•  V  » ~  *  » 

,*  v 

V'.-V.^V 

-N  >  : 

'  s-  s-  s' 

•  » 


Figure  22.  Crack  patterns  of  models  Cl  and  C2 


expected  to  occur  at  one  of  the  crown  sections.  The  formation  of  the  first 
hinge  did  occur  at  the  crown  such  that  the  tension  concrete  was  cracked  and 
the  compression  concrete  crushed.  Then  plastic  hinges  were  formed  at  the  two 
springing  line  sections  formir  ,  a  three-pinned  arch  structure.  The  last  hinge 
was  formed  following  a  secondary-type  compression  failure  at  the  other  crown 
section  forming  the  four-hinge  collapse  mechanism.  The  final  crack  patterns 
are  shown  in  Figure  23  with  the  observed  order  of  cracking  and  crushing  of 
concrete  labelled. 
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Acoustic  Emissions 


40.  After  the  test  on  Cl  was  completed,  it  was  suggested  that  it  may 
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measured  as  accelerations  occurring  above  a  threshold  acceleration  value) 
versus  load.  In  Figure  24,  the  location  of  the  transducer  is  superimposed  on 
the  crack  pattern  for  models  C2  and  C3.  In  Figure  25,  the  number  of  counts 
are  plotted  for  the  first  and  second  (retest)  test  of  C2.  Note  that  the 
vertical  scale  for  the  number  of  counts  is  a  factor  of  10  smaller  for  the 
retest  than  the  first  test.  In  Figure  26,  the  count-load  plot  is  shown  for 
the  test  of  C3  to  failure.  Distinct  differences  in  the  count-load  plots  are 
observed  between  a  failed  specimen  (C3)  and  an  intact  specimen  (C2). 
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Figure  26.  Count-load  plot  for  test  of  model  C-3  to  failure 


Load  Distributions 


42.  Data  from  the  load  cells  are  presented  in  radial  load  distribution 
plots.  Load  distributions  at  a  given  time  are  plotted  on  the  same  axis  to 
give  the  picture  of  the  growing  load  distribution  with  increasing  pressure. 
Star  symbols  represent  the  predicted  load  distribution  for  the  last  load  step 
Figures  27,  28,  and  29  contain  examples  of  the  load  distributions  recorded. 
All  the  loads  check  to  within  0.5  to  4  percent  of  static  equilibrium.  Equi¬ 
librium  checks  were  performed  on  the  measured  load  distributions  to  evaluate 
errors  in  computing  static  equilibrium.  Looking  at  Figure  27  (load  distribu¬ 
tions)  the  equilibrium  checks  were  performed  by  summing  all  negative  vertical 
and  horizontal  force  components  and  comparing  the  sums  to  all  positive  ver¬ 
tical  and  horizontal  force  components.  The  differences  range  from  0.09  to 
4.56  percent,  as  shown  in  Table  6. 
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Figure  29.  Load  distributions  recorded, 
specimen  C3 


Load-Strain  Plots 


43.  The  data  from  the  strain  gages  on  the  steel  reinforcement  and  con¬ 
crete  are  plotted  as  load  versus  strain.  Averaged  load-strain  plots  for 
tests  of  C2  and  C3  are  shown  in  Figures  30  through  33.  The  load-strain  plots 
for  each  good  gage  channel  for  tests  of  Cl,  C2,  and  C3  are  shown  in  Appen¬ 
dix  B.  As  mentioned  earlier,  the  positive-  and  negative-moment  sections  refer 
to  the  crown  and  springing  line  sections,  respectively. 

Strain  Distributions 


44.  Strains  were  plotted  with  all  gages  in  the  two  positive-moment 
sections  lumped  into  one  section  (assuming  symmetry  of  loading  and  model  de¬ 
formation).  Likewise,  for  the  negative-moment  and  the  zero-moment  sections 
(the  45-deg  sections),  the  strains  were  plotted  in  one  section.  This  gives  a 
graphical  view  of  the  scatter  in  the  data. 

45.  The  strain  distributions  of  the  positive-  and  negative-moment  sec¬ 
tions  were  fitted  with  both  a  hyperbolic  and  a  linear  curve.  The  formulations 
for  the  curve  fits  are  included  in  Appendix  C.  A  least-squares  criterion  was 
used  for  the  curve  fit.  Figures  34  and  35  are  examples  of  the  curve  fits 
representing  the  strain  distributions  in  the  positive-  and  negative-moment 
sections,  respectively,  in  C2  for  the  same  time  load  step  (or  time  increment). 
Figure  36  is  a  typical  plot  of  the  strain  distribution  at  the  45-deg  section 
(theoretically  the  zero-moment  section). 

46.  In  Table  7  the  results  from  the  curve  fits  are  listed.  The  re- 

p 

suits  in  the  column  labelled  r  are  the  values  of  the  correlation  coef¬ 
ficient  squared.  This  value  is  a  measure  of  the  scatter  in  data  from  the  as¬ 
sumed  strain  distribution  (linear  or  hyperbolic).  The  two  values  solved  by 
the  curve  fits  are  the  location  of  the  neutral  axis,  c  ,  or  the  depth  of  the 
compressive  zone  and  the  maximum  compressive  strain,  eu  . 

Resolving  Moments  and  Thrusts  from  Data  Fit 

47.  Moments  and  thrusts  were  computed  for  each  strain  distribution 
using  the  section  designs  in  Figures  17  and  18  and  the  material  properties  in 
Table  4.  These  values  were  used  as  input  to  a  computer  program  of  the  curved 
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Figure  33.  Average  of  load-strain  plots: 
(1)  gages  9  and  10;  (2)  gages  11  and  12; 
(3)  gages  13  and  16 
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34.  Strain  distribution  in  positive-moment  section 
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Figure  36.  Typical  strain  distribution  at  45-deg  sections 
(zero-moment  sections) 

beam  analysis.  Moments  and  thrusts  were  computed  from  both  the  hyperbolic  and 
linear  strain  distributions  once  the  neutral  axis  location  and  the  maximum 
compressive  strain  were  known  with  the  input  of  the  material  properties.  For 
the  curved  beam  analysis  outlined  in  Part  II,  either  strain  distribution  may 
be  assumed  to  compute  the  moments  and  thrusts.  In  Table  8,  all  the  moments 
and  thrusts  are  listed  for  different  load  levels.  For  each  load  level,  three 
comparisons  for  moments  and  thrusts  are  made  among:  (a)  curved  beam  theory 
(application  to  RC  curved  beam),  (b)  straight  beam  theory,  and  (c)  a  1-D 
linear  elastic  structural  analysis. 

48.  Values  from  Table  8  are  plotted  on  the  interaction  diagrams  (plot 
of  thrust  and  moment  capacity  of  a  section)  shown  in  Figures  37  to  41.  The 
five-point  loading  load  path  is  the  predicted  1-D  elastic  response  due  to 
five-point  loads  per  quarter  of  the  model.  The  load  path  marked  by  p(e)  is 
the  predicted  1-D  elastic  response  due  to  a  continuous  normal  pressure.  The 
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Figure  37.  Resolved  moments,  thrusts,  and  load  paths, 
specimen  Cl,  positive-moment  section 


•  ■ 

•v 


M 

Sl-vV 

■».  ■>-.  • .  - 
-•.  -w  v 
-•.v.v.  - 

V.  V 

v.v>  .v 

V.V.V^ 


•  _ 

T,  •  ,*  »  *1 

•  >;-v 


>» 

\  % 

V."  v 
A.' 


■*#*.'*1 

■>>.v 

■vvv 


.".V 


Hi 


•  -•  .•  •  •  -  •  .  .  •  .->•  •  •  . 
Awlv’  •.  .  „  • v‘.  v'i'  V  '-  .  • v’_v .  O  v*.  v'_v' 


1.0 


M 

fc'  bhz 


Figure  38.  Resolved  moments,  thrusts,  and  load  paths, 
specimen  C2,  positive-moment  section 
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Figure  40.  Resolved  moments,  thrusts,  and  load  paths 
specimen  C3,  positive-moment  section 
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Figure  41.  Resolved  moments,  thrusts,  and  load  paths, 
specimen  C3,  negative-moment  section 

open  circular  and  triangular  symbols  marked  with  numbers  near  the  predicted 
load  paths  represent  resolved  moments  and  thrusts  for  a  particular  reference 
load  level  as  computed  by  assuming  a  straight  or  curved  beam  response.  The 
interaction  diagrams  have  been  computed  assuming  a  curved  beam  response.  In 
the  interaction  diagrams,  the  triangle  symbols  connected  with  solid  lines 
represent  the  ultimate  curve  and  the  solid  symbols  represent  the  yield  moment 
curve. 

49.  Data  for  the  negative-moment  section  of  Cl  were  recorded  on  another 
analog  32-track  tape  and  were  not  complete  because  of  a  malfunction. 


PART  VI:  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


50.  The  first  two  specimens,  Cl  and  C2,  did  not  fail  when  tested  to  the 
full  capacity  of  the  loading  apparatus.  Specimen  C3  failed  in  flexure  ini¬ 
tially  and  ultimately. 

51.  The  first  test  on  Cl  indicated  no  visible  signs  of  iistress  on  the 
model.  After  the  test,  it  was  discovered  that  one  of  the  hydraulic  circuits 
had  inadvertently  been  shut  off;  thus,  all  of  the  load  that  was  believed  to  be 
on  the  specimen  was  not.  The  hydraulic  circuit  for  jacks  8A  and  8B  was  closed 
(see  Figure  21).  With  the  absence  of  the  loads  in  8A  and  8B,  there  is  a 
slight  skewness  to  the  load  distribution.  The  pressure  in  this  circuit  is 
90.2*1  percent  of  the  maximum  pressure  in  circuit  1  (jacks  1A  and  IB).  At  the 
crown,  the  effect  of  the  load  in  jacks  8A  and  8B  accounts  for  8.5  percent  of 
the  total  thrust  and  37.1  percent  of  the  total  moment.  To  account  for  the 
absence  of  load  in  jacks  8A  and  8B,  the  entire  conduit  response  is  analyzed 
and  the  results  shown  in  Appendix  A  (Figure  A 2). 

52.  Specimen  C2  was  by  design  a  replica  of  Cl.  The  purpose  of  conduct¬ 
ing  two  identical  tests  was  to  establish  consistency.  Retests  on  Cl  (and  C2) 
established  consistency  within  a  specimen.  Because  of  the  slight  discrepancy 
between  the  initial  load  distributions  of  Cl  and  C2,  the  tests  were  not 
exactly  the  same.  For  the  first  test  on  C2,  all  pressure  channels  were 
checked  to  assure  that  all  loads  were  acting  on  the  model.  Specimen  C2  was 
not  failed  either  when  tested  to  full  capacity  of  the  loading  system.  After 
each  retest  (three  on  Cl,  one  on  C2)  more  cracks  were  visible. 

53.  Retests  conducted  on  Cl  and  C2  evaluated  the  strength  of  each  model 
through  a  cycle  of  loading.  Even  after  a  reloading  of  the  models,  no  signifi¬ 
cant  damage  was  observed. 

54.  Specimen  C3  was  the  only  model  loaded  to  ultimate  failure.  The 
progression  to  failure  was  described  in  the  test  observations  in  Part  IV.  A 
flexure  collapse  mechanism  formed  such  that  "plastic  hinges"  developed  at  the 
crown,  invert,  and  the  two  springing  lines.  No  signs  of  shear  failure  were 
observed.  Specimen  C3  was  the  thinner  model  with  R/h  =4.0  . 

55.  In  Table  7,  all  the  results  of  the  fits  to  the  strain  distributions 

p 

are  listed.  The  values  of  r  (the  correlation  coefficient  squared)  ranged 
from  0.749  to  0.989  for  the  hyperbolic  fit  and  0.666  to  0.985  for  the  linear 

p 

fit.  The  value  of  r  implies  what  percentage  of  the  scatter  in  the  data  can 


be  represented  by  the  specific  curve  fit.  The  results  in  Table  7  show  no  sig¬ 
nificant  difference  between  the  hyperbolic  and  the  linear  fit.  However,  these 
results  are  based  on  data  from  three  specimens,  which  is  a  limited  statistical 
sample.  Most  of  the  strain  gages  were  placed  in  two  locations  and  some  in 
three.  More  meaningful  strain  gradients  should  be  obtained  to  better  quantify 
the  statistical  comparison  between  the  hyperbolic  and  linear  fit.  This  could 
be  accomplished  by  placing  more  circumferential  strain  gages  across  the  sec¬ 
tion  on  future  tests.  Accordingly,  more  tests  results  will  be  needed  to  con¬ 
clude  whether  the  current  straight  beam  analysis  is  appropriate  for  conduits 
with  R/h  as  low  as  2.5. 

56.  Table  8  and  Figures  37  to  Ml  contain  the  results  of  the  resolved 
moments  and  thrust  for  specimens  Cl,  C2,  and  C3.  In  Figure  37,  the  moment s 
and  thrusts  form  a  straight  load  path  and  are  contained  within  the  interaction 
diagram.  In  Table  8,  the  experimental  values  are  compared  to  theoretical 
values  (computed  from  the  1-D  structural  analysis).  For  specimen  Cl, 
differences  in  thrust  range  from  36  to  64  percent  between  the  linear  fit  and 
the  predicted  value  (theoretical).  The  differences  in  the  moments  ranged  from 
46  to  61  percent.  The  largest  differences  occurred  at  the  lowest  load 
levels.  These  results  are  for  the  internal  forces  resolved  at  the  crown. 

57.  For  specimen  C2  at  the  crown,  differences  in  the  thrust  values 
ranged  from  0.51  to  26  percent,  with  the  moment  values  differing  from  30  to 
40  percent.  At  the  springing  line,  thrusts  differed  from  7.9  to  35  percent 
and  moments  differed  from  24  to  61  percent.  In  Figure  38,  the  experimental 
values  plot  a  consistently  different  load  path  than  predicted.  The  moment  and 
thrust  at  the  last  load  level  plots  are  within  the  interaction  diagram  and 
this  is  consistent  with  the  observed  experimental  behavior.  Specimen  C2  did 
not  fail  as  predicted  by  the  five-point  loading  path  (the  predicted  loading 
path  is  theoretically  between  the  p(8)  and  five-point  loading  paths).  Even 
after  the  observed  first  crack  the  load  path  appears  to  be  undisturbed.  In 
Figure  39,  the  plotted  experimental  values  of  moment  and  thrust  appear  to 
agree  closely,  with  the  predicted  values.  Even  though  the  experimental  and 
theoretical  load  paths  agree  closely,  there  is  a  difference  between  the  values 
at  specific  load  levels.  No  signs  of  failure  existed  on  specimen  C2  at  the 
springing  line,  which  is  consistent  with  plotted  results. 

58.  For  specimen  C3,  the  only  specimen  tested  to  failure,  the  thrusts 
and  moments  at  the  crown  differed  from  12  to  35  percent  and  23  to  57  percent, 


respectively.  At  the  springing  line,  the  thrusts  differed  from  14  to 
24  percent  and  the  moments  differed  from  33  to  54  percent.  In  Figure  40,  the 
experimental  and  theoretical  load  paths  at  the  crown  are  plotted.  Again  the 
difference  between  the  two  load  paths  is  observed.  The  structure  yielded  at 
about  14.5  kips  (maximum  load  in  jack  1A)  as  seen  in  the  load-strain  plots  in 
Appendix  B  (strain  gages  3  and  4)  and  as  shown  by  the  intersection  of  the  ex¬ 
perimental  load  path  with  the  yield  line  of  the  interaction  diagram.  The  pre¬ 
dicted  ultimate  and  yield  moments  are  very  close,  as  shown  in  Figures  40  and 
41.  However,  a  closer  look  at  the  load-strain  plots  on  pages  B55  and  B56  (of 
Appendix  B)  reveal  that  the  ultimate  load  resistance  occurred  at  approximately 
the  32,500-lb  reference  load  level.  This  is  about  twice  the  first  yield  load 
level.  Reference  to  the  load-strain  plots  of  the  tension  steel  at  the  crown 
(pages  B45,  B46,  B49,  and  B57)  indicates  that  the  first  yield  occurs  at  about 
the  15,000-lb  reference  load  level.  Second  yield  occurs  at  the  springing  line 
(pages  B51,  B52,  and  B59).  The  plots  on  pages  B55  and  B56  are  the  load-strain 
results  of  the  concrete  compressive  strains  at  the  intrados  of  the  springing 
line.  The  second  and  third  plastic  hinges  formed  to  complete  a  three-pinned 
arch  structure  before  final  collapse  due  to  the  fourth  hinge. 

59.  In  Figure  41  it  is  shown  that  the  capacity  of  the  springing  line 
section  was  not  exceeded  at  15  kips.  This  was  observed.  The  analytical  and 
experimental  load  paths  at  the  springing  line  section  agree  very  closely,  but 
the  values  at  specific  load  levels  differ  as  discussed  previously. 

60.  Three  consistent  observations  among  the  tests  are:  (a)  the 
discrepancy  between  analytical  and  experimental  load  paths  at  the  crown,  (b) 
agreement  between  analytical  and  experimental  load  paths  at  the  springing 
lines,  and  (c)  comparison  of  observed  cracking  failures  with  predicted 
location  of  cracking. 

61.  Three  effects  will  be  discussed  that  can  account  for  the  observed 
differences  between  predicted  and  experimentally  computed  internal  forces  at 
the  crown  and  the  springing  line  sections.  First,  the  tensile  capacity  of  the 
concrete  has  been  ignored.  Looking  at  Figures  42  and  43  (reprints  of  Fig¬ 
ures  34  and  35,  the  strain  distributions  at  10  kips)  the  maximum  thrust  and 
moment  due  to  tension  carried  by  concrete  are  computed  using  the  linear  fit 
for  convenience.  The  moments  are  seen  to  be  affected  more  than  the  thrusts. 
However,  in  this  test,  10  kips  is  a  low  load  level  and  the  effect  of  the  ten¬ 
sile  capacity  of  the  concrete  would  be  expected  to  be  less  for  the  higher  load 


"A  biaxial  stress  condition  occurs  if  the  principal 
stresses  act  only  in  two  directions;  that  is,  the 
stresses  act  in  one  plane  and  the  third  principal 
stress  is  zero. .. [Figure  42]  presents  the  combina¬ 
tion  of  direct  stress  in  two  directions  which  caused 
failure  of  concrete  subjected  to  biaxial  compression 
may  be  as  much  as  27  percent  higher  than  the  uniaxial 
strength.  For  equal  biaxial  compressive  stresses, 
the  strength  increase  is  approximately  1 6  percent. 

The  strength  under  biaxial  tension  is  approximately 
equal  to  the  uniaxial  tensile  strength.  Note,  how¬ 
ever,  that  combined  tension  and  compression  loadings 
reduce  both  the  tensile  and  the  compressive  stresses 
at  failure"  (from  Park  and  Paulay  (1975)). 


Figure  44.  Biaxial  strength  of  con¬ 
crete,  f  =  uniaxial  strength 
(Park  and  Paulay  1975) 


when  the  load-strain  plots  are  compared  betweeen  the  initial  and  the  reloading 
test  of  C2  in  Figures  45-48.  In  all  the  initial  load-strain  plots,  a  non¬ 
linear  response  is  observed.  In  all  the  load-strain  plots  of  the  corres¬ 
ponding  reloaded  channels,  it  appears  that  the  structure  is  behaving  linearly 
elastic.  The  structure  was  not  reloaded  past  the  maximum  load  of  the  initial 
test.  According  to  the  theory  of  plasticity,  after  an  initial  loading  of  a 
structure  (or  material)  past  yield,  the  structure  will  behave  in  an  elastic 
manner  until  loaded  past  yield  again  (Anderson,  Haelsig,  and  Reifel  1966). 
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07/14>82  R0978 


b.  RELOADING 


C2.2  61  7S2 
S-1 
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CHANNEL  NO.  21  20167  1 

07/16/82  R0094 
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Figure  45.  Comparison  of  load-strain  plots, 
initial  test  and  reloading 
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b.  RELOADING 
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Figure  46.  Comparison  of  load-strain  plots, 
initial  test  and  reloading 
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PART  VII:  FUTURE  WORK 

64.  In  FY  83,  model  RC  circular  conduits  were  planned  for  testing  to 
ultimate  failure.  The  results  of  these  tests  will  provide  qualitative  infor¬ 
mation  about  failure  modes  and  will  be  the  subject  of  a  future  report. 

65.  Finite  element  studies  will  also  continue,  using  both  linear  and 
nonlinear  models.  The  effect  of  the  biaxial  response  of  concrete  will  also  be 
studied. 

66.  The  problem  of  shear  capacity  of  the  RC  conduits  will  be  ad¬ 
dressed.  New  criteria  will  be  evaluated  and  revised  if  necessary. 
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PART  VIII:  CONCLUSIONS 


67.  Initial  curvatures  of  Corps  curvilinear  conduits  do  affect  their 
flexural-axial  resistance,  but  when  compared  to  straight  beam  behavior,  the 
effect  of  the  initial  curvature  is  apparently  not  significant.  It  has  also 
been  shown  from  curve  fits  to  experimental  strain  data  that  there  is  statis¬ 
tically  no  difference  between  hyperbolic  and  linear  strain  distributions  in 
circular  conduits.  More  tests  will  be  needed  to  obtain  a  better  statistical 
sample  and  investigate  even  more  extreme  curvatures  than  studied  in  this 
report.  From  future  tests  it  will  be  determined  whether  the  normal  design 
practice  for  computing  the  moment-thrust  resistance  of  uncurved  members 
provides  an  appropriate  safety  margin  for  Corps  cross  sections.  The  shear 
resistance  of  these  Corps  structures  was  not  considered  herein;  therefore,  it 
is  the  subject  of  a  further  investigation. 

68.  A  nonlinear  finite  element  model  for  concrete  was  shown  appli¬ 
cable.  Also  a  1-D  structural  analysis  and  2-D  Airy  stress  solutions  compared 
favorably  for  ev.-uuating  the  response  of  the  elastic  behavior  of  RC  conduits 
under  static  pressure  loading.  Three  RC  conduit  specimens  were  loaded  and 
these  tests  established  a  standard  testing  procedure.  Of  the  three  specimens 
only  one  failed  ultimately.  Cracking  occurred  during  tests  on  specimens  Cl 
and  C2.  The  third  specimen  formed  a  flexural  collapse  mechanism  as  was 
predicted.  The  resultant  moments  and  thrusts  from  the  data  differ  from  the 
internal  forces  predicted  from  the  linear  elastic  1-D  structural  analysis. 
Three  reasons  cited  for  these  differences  are: 

a.  The  tensile  capacity  of  concrete  was  neglected  (this  affects 
results  at  the  smaller  load  levels). 

b.  Biaxial  response  of  concrete  was  not  accounted  for  in  the 
constitutive  model. 

c.  Plastic  behavior  of  the  structure  was  not  accounted  for  in  the 
structure  analysis. 

Items  b  and  c  are  the  two  most  important  effects  to  be  studied  for  future 
analysis. 
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Table  1 

Characteristics  of  Corps  Curvilinear  Conduits 


f  ' 
c 

Pg 

f 

y 

R 

d 

d' 

e 

Structure 

Shape 

ks  i 

% 

% 

ksi 

h 

h 

h 

h 

Liu  (1980) 
Examples  1,2 

Circular 

4 

0.72 

0.36 

40 

3.0 

0.875 

0.125 

0.417 

Liu  (1980) 
Examples  3,4 

Oblong 

4 

0.72 

0.36 

40 

3.0 

0.875 

0. 125 

0.289 

Liu  (1980) 
Example  5 

Horseshoe 

4 

0.24 

0 

40 

3.4 

0.786 

-- 

0.239 

Mattock, 

Kriz,  and 
Hognestad 
(1961) 

Circular 

4 

0.95 

0.36 

40 

3.8 

0.75 

0.25 

Mattock, 

Kriz,  and 
Hognestad 
(1961) 

Circular 

4 

0.82 

0.27 

40 

3.0 

0.81 

0.19 

Hognestad, 

Hanson, 

and 

McHenry 

(1955) 

Circular 

4 

0.44 

0.22 

40 

2.5 

0.88 

0.  10 

Note:  fg  =  design  28-day  uniaxial  compressive  concrete  strength;  fy  =  design 

yield  strength  of  steel;  p  =  total  A  /bh  ,  gross  steel  ratio; 

g  s 

o'  =  A'/bh  ,  gross  compression  steel  ratio;  h  =  depth  of  section; 
go 

R  =  radius  to  middepth  of  section;  e  =  eccentricity  =  M/P  (moment/ 
thrust);  d  =  effective  depth  of  tension  steel;  d'  =  depth  to 
compression  steel. 


Hydraulic  System 
Hydraulic  cylinder  jacks 

Hydraulic  power  supply 

Hydraulic  controller 
(pressure  regulator) 

Instrumentation 
Load  cells 

Strain  gages 
Deflection  gages 


Double  Acting 

Hydroline  Series 

Hydroline  Cylinders 

Model:  N2R,  Serial  Suffix:  01R, 

Bore — 5  in.,  Stroke — 6  in. 

Rated  for  58.9  kips  at  3,000  psi 

Built  by  Activation,  Inc.,  Model  #1-1671-1 
Serial  No.  80-1751  rated  for  3,000  psi 

Hydraulic  Load  Maintainer 

Challenger  Model  10M 

Serial  No.  5M1142  rated  for  5,000  psi 


Tnterface  Model:  1220  AR50K,  Precision 
Universal  Load  Cells,  50  kips 

Micro  Measurements 

Linear  Variable  Differential  Transformer 
(LVDT) 


2  tape  recorders 


32-track  Sabre  IV 


Table  5 

Load  Verification  Results 
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Table  6 

Load  Equilibrium  Checks 
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P  sin  0  sin  $  ,  41  <  0 


P  jsLn  8  sin  $  -  cos  (41  -  9)J  ,  <t>  _> 


where 

N($)  =  internal  thrust 
N ( 4> )  =  internal  moment 
V(< |>)  =  internal  shear 
Mc  =  crown  moment 
Mg  =  springing  line  moment 
Rc  =  crown  thrust 
Rg  =  springing  line  thrust 
<t>  =  angle  from  springing  line  to  crown 
0  =  angle  at  which  P  is  applied 
P  =  normal  load  applied  at  0 
r  =  radius  of  1-D  structure 

4.  The  resultant  internal  forces  as  given  by  a  previous  solution 
(Anderson,  Haelsig,  and  Reifel  1966)*  to  the  structure  in  Figure  15  (in  the 
main  text)  under  the  continuous  pressure  distribution  are  given  by 


wn  k  .  v 
M(0)  =  -g2-  (Kv  -  Kj  cos  29 

"<«>  =  5  »0  [(2KL  *  Kv)  *  (Kv  -  kl) 


2 

COS  0 


where 


v(e)  =  i  (kl  - 


sin  29 


0  =  angle  measured  from  the  springing  line 
w  =  maximum  pressure  (pmax) 


Rq  =  outside  radius 


R  =  middepth  section  radius 
=  lateral  pressure  coefficient 
Kw  =  vertical  pressure  coefficient 


he  References  at 


1.5 


for  3:1  loading  =>  KL/Ky  =  1/3  =>  KL  =  0.5  ,  Ky  = 

5.  The  EM  1110-2-2902  (Headquarters,  Department  of  the  Army  1969) 
design  loadings  give 


and 


M(e) 
N(e) 
V(  9 ) 


wR  R 

=  -f-  (Kv  -  kl)  cos  26 


2  2 
wR  /  sin  0  +  Ky  cos  9 


•) 


wR 


(kl  -  Kv) 


sin  28 


6.  The  resulting  internal  forces  from  the  Airy  stress  solution  to  the 
structure  in  Figure  15  (in  the  main  text)  under  the  continuous  pressure 
distribution  are 


M(8)  = 


where 

e  =  r  -  — 

j  jft  (b2  -  a2)  ♦  3B  (b1*  -  a1*)  -  3C 

- cos  28  +  F  ■  H[tbtb)  -  in(aj]j 

ra 

jpA(b  -  a)  *  4B  (b3  -  a3) 

'  2C  (bT-^)]  -  a)  .  hQ  -  I)J 

where 


b  =  outer  radius 


a  =  inner  radius 


B  = 

'“I 

2  2 
(d  ♦  3a  ; 

C  = 

p  ,,  U  6,2, 

a  q  (3a  b  +  a  b  ) 

F  = 

-3a 

2 

.2  2  c 
b  -  a 

H  r 

6aa2 

2 

,2  2  * 
b  -  a 

D  = 

,  4  2  .2  4  6, 

(2b  a  +  b  a  +  a  ) 

‘m  = 

R 

p  =  varying  pressure  component  magnitude, 


(K  +  Ky)w 

q  =  uniform  pressure  component  magnitude  - r- — — 


7.  Accounting  for  the  absence  of  the  load  at  locations  8A  and  88,  the 
effect  of  this  absence  is  solved  by  looking  at  the  entire  conduit  (1-D)  loaded 
by  a  two-point  load  and  subtracting  the  effects  contributed  at  the  crown  and 
the  springing  line.  In  Figure  A2,  the  results  of  the  classical  two-point 
loading  on  a  ring  (Timoshenko  1941)  are  shown. 
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APPENDIX  C 

STRAIN  DISTRIBUTIONS 
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Nonlinear  Least-Squares  Regression  Analysis 
(Hyperbolic  Strain  Distribution) 


1.  The  hyperbolic  strain  distribution  according  to  Timoshenko's  (1941) 
curved  beam  theory  is  (for  the  positive  moment  case): 


,  .  y  (P  +  l)e2 

6  y  ~  (a  h  \ 

[*  +  2  -  c  *  y) 


Discretely,  replace  e(y)  with  e^y^)  or  and  y  with  y ^  . 
Therefore  the  data  will  be  curve  fitted  by  the  equation: 


R'y.e 
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i  "  (R*  -  c  +  yi)c 


where 


R'  =  (R  +  §) 

e  and  c  =  undetermined  coefficients  to  be  solved  for  by  the  least- 
u  squares  method 

2.  The  error  term  ei  is: 
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ei  "  ei  “  (R'  -  c  +  y^c 


The  sum  of  the  squares  of  the  error  term  is: 
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where  N  is  the  number  of  data  pairs  (ei,yi) . . .(e^, yN ) .  Let  eu  =  a  * 
maximum  compressive  strain  and  c  =  B  ,  location  of  the  neutral  axis, 
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6.  Therefore,  using  Equations  C6  and  C7,  the  following  method  was  used 

to  solve  for  a  and  6  : 

a.  Choose  a  starting  value  for  8  (a  realistic  value  is  chosen 
using  an  "eyeball"  fit  to  the  data  first). 

b.  Calculate  a  from  Equation  C6  where  6  is  computed  from  a. 

c.  Rewrite  Equation  C7  as: 
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d.  Using  a  root  finder,  determine  a  new  value  of  B  and  repeat 

steps  a  through  d  until  a  convergence  criterion  is  met  such  that 


g.,(8)  -  g2(8) 


<  t 


where  g1  and  g2  are  two  function  values  from  successive  iterations  solving 
for  a  and  8  and  t  is  a  user-specified  tolerance.  There  are  practical 
limitations  for  the  value  of  t  ;  the  limit  of  t  becomes  clear  when  the 
above  solution  scheme  is  programmed  and  run  on  the  computer.  The  solution  for 
a  and  B  for  obtaining  the  hyperbolic  curve  fits  to  the  strain  data  was 
obtained  on  a  4081  Tektronix  computer. 

7.  Similarly,  using  the  aforementioned  procedure  by  substituting  for 
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8.  Equations  C9  and  CIO  were  derived  from  Equations  C11  and  C12 
below.  The  hyperbolic  strain  distribution  for  a  negative  moment  section  is: 
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and  where  R'  =  R  -  h/2  ,  the  sum  of  the  squares  of  the  errors  (or  deviations) 
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and  letting  a  =  eu  and  B  =  c  , 
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Correlation  Coefficient 
(Miller  and  Freund  1977) 

9.  The  correlation  coefficient  squared  gives  statistical  inference  to 
quantify  the  correlation  between  the  data  and  the  fitted  curve.  The  equation 
below  measures  the  correlation  by: 
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where 

p 

r  =  the  correlation  coefficient  squared 
e :  =  measured  value 
e |  =  predicted  value 
e  ^  =  average  of  measured  values 

10.  The  data  that  were  fitted  with  the  hyperbolic  curve  were  also 
analyzed  by  the  linear  regression  method.  This  gave  a  comparison  of  the 
statistical  differences  between  a  hyperbolic  and  a  linear  strain  distribution 
for  the  same  data.  The  linear  regression  formulas  can  be  found  in  many  sta¬ 
tistics  references.  Some  hand  calculators  also  have  linear  regression 
equations  preprogrammed  by  keystroke  operations. 
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APPENDIX  D 


DERIVATION  OF  K-,  ,  Kz  ,  K3  CURVES  IN  FIGURE  8 
DEPENDENT  ON  CYLINDER  STRENGTH 
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Assuming  eQ  =  0.002  and  eu  =  0.00H  - 
the  following  is  derived: 


,  where  f’  is  ksi  units  . 
c 


Positive  Moment  Section.  Neutral  Axis 
at  Inner  Radius  (r  =  R  -  h/2) 
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,  where  Z  =  h  +  y 
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Negative  Moment  Section.  Neutral  Axis  at 
Outer  Radius  (r  =  R  +  h/2) . 


(MU) 

(8-M) 


Derivation  of  K1  ,  K2  ,  {refer  to  Figure  D1) 


K.  is  the  ratio  of  average  stress  to  maximum  stress 
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K2  is  the  ratio  of  depth  to  resultant  at  concrete  and  depth  to  neutral  axis 
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for  positive  moment  section 


'A.vTy.v.Mv.vv  v.wv  v  h.*  o  u  W 

>■-  • . -vy.  ■f.  ■ 

w.\  /,>■, 

. « .VA.V' 


K„  = 


E 

i-1 

E  fi 

i=1 


n 

2  IK 

i- 1  n  1 

=  1  _  -  for  negative  moment  section 

n  r 

E  fi 

i=1 

is  the  ratio  of  maximum  stress  to  cylinder  strength,  f'  (psi) 
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Inner  radius 

Width  of  section:  outer  radius  as  defined  in  Appendix  A 
Depth  of  section  in  compression 
Resultant  compressive  concrete  force 
Resultant  force  in  the  compressive  steel 

Effective  depth  of  tension  steel 
Depth  to  compression  steel 

Small  angle  between  two  neigboring  cross  sections 
Eccentricity  (M/P) 

Initial  elastic  concrete  modulus 

Modulus  of  elasticity  of  steel  (29,000,000  psi) 

Concrete  stress 

28-day  compressive  strength  of  standard  6-  by  12-in.  concrete 
cylinder 

Compressive  strength  of  concrete  in  reinforced  concrete  members 
(0.85  f£) 

Steel  stress 

Uniaxial  strength 

Yield  strength  of  reinforcing  steel 

Tensile  yield  strain 

Overall  depth  of  section 

Lateral  pressure  coefficient 

Vertical  pressure  coefficient 

Bending  moment 

Crown  moment 

Springing  line  moment 

Internal  moment 

Internal  thrust 

Varying  pressure  component  magnitude 
Centric  thrust 

Uniform  pressure  component  magnitude 

Radius  of  initial  curvature  to  neutral  axis  of  section 

Radius  of  initial  curvature  to  middepth  of  section 


Crown  thrust 
Outside  radius 


Springing  line  thrust 

Resultant  force  in  the  tensile  steel 

Internal  shear 

Maximum  pressure 

Distance  from  the  neutral  axis 

Small  angle  of  rotation  due  to  moment  and  thrust 

Strain 

Maximum  concrete  strain 

Strain  at  which  the  maximum  concrete  stress  f”  occurs  (0.002) 

Tensile  reinforcement  strain 

Limiting  useful  concrete  strain  (0.003) 

Angle  at  which  P  is  applied,  or  angle  from  springing  line  of 
section  of  interest 

Gross  tension  steel  ratio 

Gross  compression  steel  ratio 

Angle  from  springing  line  to  section  of  interest 


